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By considering the constrained motion of classical spins in a geometrically frustrated magnet, we 
find a dynamical freezing temperature below which the system gets trapped in metastable states 
with a "frozen" moment and dynamical heterogeneities. The residual collective degrees of freedom 
are strongly correlated, and by spontaneously forming aggregates, they are unable to reorganize 
the system. The phase space is then fragmented in a macroscopic number of disconnected sectors 
(broken ergodicity), resulting in self-induced disorder and "thermodynamic" anomalies, measured by 
the loss of a finite configurational entropy. We discuss these results in the view of experimental results 
on the kagome compounds, SrCr9pGai2_9pOi9, (H30)Fe3(SO4)2(OH)6, Cu3V207(OH)2.2H20 and 
Cu3BaV208(OH)2. 



I. INTRODUCTION 

Certain magnetic compounds lack conventional mag- 
netic long-range order but develop static order below a 
temperature Tg, with locally "frozen" spins. Well-known 
examples are spin-glasses but there are now examples of 
geometrically frustrated compounds with somewhat dif- 
ferent microscopic properties. They are dense (one spin 
per site on a periodic lattice), but have a rather small 
"frozen" moment. 

The physical origin of such glassy-like phases is an in- 
teresting issue. It may be a spin-glass phase associated 
with weak quenched disorder—"— or it may be more intrin- 
sic to the pure compound and its geometrical frustration. 
For example, structural glasses do lack quenched disorder 
but are out-of-equilibrium with a relaxation time longer 
than the observation time of the experiment. It is a gen- 
eral idea that the frustration, by suppressing long-range 
order, may lead to glassy-like phases.—"^ In the present 
paper, we study the relaxation to equilibrium of the dy- 
namics of a simple spin model, in the presence of geo- 
metrical frustration. We find that the spin relaxation 
is nonexponential and develops spontaneously two time- 
scales below a crossover temperature Td- The system 
does not slow down uniformly in space; instead, it devel- 
ops some fast-moving and slow-moving regions character- 
ized by an emergent length-scale (called "dynamical het- 
erogeneities" in the context of structural glasses^) . Below 
a second crossover temperature Tg, the slow- moving spins 
may appear "frozen" on the experimental time-scale, i.e. 
the system has fallen out-of-cquilibrium. In this case, the 
system is found to be trapped into one of an exponential 
number of metastable states and some local disorder is 
self-induced. 

Competing local spin interactions resulting e.g. from 
the geometrical frustration of the lattice tend to sup- 
press the magnetic long-range order. At low tempera- 
tures, some local correlations appear and the system is 
in a collective paramagnetic regime. The spin dynamics 
is different from that of a high-temperature paramagnet: 
the system still has a macroscopic number of accessible 
states but these states are locally constrained. The spin 



dynamics is hindered by these local constraints: single 
spin flips become suppressed if they violate local arrange- 
ments and the degrees of freedom acquire a more col- 
lective nature, which, in the present context, are loops 
(or "strings") of spins. The issue is whether these co- 
operative excitations are efficient enough to reorganize 
the system as in the liquid state (here the paramagnetic 
state) or if the system is "jammed" . Such excitations 
are rather ubiquitous and appear in different contexts, 
e.g. ice and fcrroelectricsj^ Stringlike excitations have 
been also identified in molecular dynamics simulations 
of structural glasses^ and were argued to indeed play a 
role in the glass transition problem)^ Here we study how 
these excitations self-organize in a simple degenerate spin 
model on a lattice, and how they permit or not, depend- 
ing on temperature, the relaxation to equilibrium. We 
find that, while the motion of long loops is very efficient 
at high temperatures, it is too slow at low temperatures 
and the residual "rapid" degrees of freedom do not lead 
to thermodynamic equilibrium. 

The magnetic materials we have in mind are highly 
frustrated systems with spins on the sites of the two- 
dimensional kagome lattice, but some spin-ice sys- 
tems on the three-dimensional pyrochlore lattice have 
a rather similar phenomenology^iiS and sustain sim- 
ilar loop excitations^^ The kagome systems have a 
spin freezing transition at Tg but the "frozen" moment 
is rather small and the system retains some dynam- 
ics below Tg. This is the case of the rather dense 
kagome bilayer SrCrgpGai2-9pOi9 (SCGO)r^ which was 
argued originally to be an unconventional spin-glass 
because (i) the specific heat is in T^jHii^ (ii) Tg is 
weakly sensitive to the chemical content p )^^i^^ (iii) 
the "frozen" moment is small and most of the sys- 
tem remains dynamicalJ^"— In the kagome hydronium 
jarosite^ (H30)Fc3(SO4)2(OH)6, the Tg does not de- 
pend much on the Fe coverage and compounds with 
100% of Fe (as the chemical formula suggests), were 
synthesized.^- Chemical disorder is certainly not ab- 
sent, though, with possible proton disorder.— Nonethe- 
less, temperature cycles below Tg, were qualitatively dif- 
ferent from that of conventional spin-glasses, and may 



point to a different nature of the phase transitioni^i^ 
More recently, two other kagome compounds were 
found: the volborthite^^ [Cu3V207(OH)2.2H20] and the 
vesignieite^l [Cu3BaV208(OH)2] . Both have a spin freez- 
ing transitioi 
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with small frozen moments j 
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NMR revealed a heterogeneous state below Tg : the NMR 
relaxation time appears to depend on the nucleus in 
volborthite, with "slow" and "fast" sites found in the 
lineshapC] ^^'^^ In vesignieite, a partial "loss" of some nu- 
clei (partial "wipeout" of the intensity) is also possibly 
indicative of sites with slower magnetic environmentsi^ 
These experiments may suggest the presence of dynam- 
ical heterogeneities 1^ These are two-dimensional sys- 
tems but a freezing transition also occurs in the hyper- 
kagome gadolinium gallium garnet, Gd3Ga50i2, a three- 
dimensional version of the kagome latticci^ However, 
not all kagome antifcrromagncts have a spin freezing 
transition. Some have antiferromagnetic long-range or- 
der, such as those of the jarosite famil}*^^ (other than 
the hydronium jarosite) or the oxalatesi^ Others may 
be quantum spin liquids, such as the herbertsmithite 
ZnCu3(OH)6Cl2 that has no phase transitions^ and a dy- 
namics down to the lowest temperatures with no clear 
energy scale in neutron inelastic scattering)^"— Such 
a broad response has some similarities with that of 
SCGOi2r— or the hydronium jarosite above the freez- 
ing temperature4S This points to competitions between 
different states and while it is possible to model some an- 
tiferromagnetic phases by appropriate interactions, e.g. 
further-neighbor interactions)^ or Dzyaloshinskii-Moriya 
interactions )^2iiS the issue of the spin freezing is delicate. 

Many theoretical studies of spin freezing phenomena 
in the context of the kagome antiferromagnet have been 
undertaken, mainly from classical or semi-classical ap- 
proaches. The role of the local collective degrees of free- 
dom (also called "weathcrvanc" modes) was put forward, 
leading to the conjecture of a spin freezing for the Heiscn- 
berg kagome antiferromagnet . ^"^'^^ It was later argued 
that distortions may help in stabilizing a "frozen" state, 
e.g. a trimerized kagome antiferromagnet has slow dy- 
namics on time-scales of single spin-flipS)^ (short com- 
pared with the time-scales probed in the present study, 
as we shall see) or distorted kagome lattices Ji It is in 
discrete spin models that a "jamming" transition was 
found, in the presence of additional interactions that fa- 
vor an ordered state: the dynamics becomes very slow 
as a consequence of a special coarsening of the domains 
of the ordered phase 4^i^ Here we shall consider similar 
discrete spins, with a different classical dynamics (not 
induced by additional interactions -the equilibrium state 
remains paramagnetic), but resulting from activated mo- 
tion within discrete degenerate states. 

The paper is organized as follows: in section |TT1 we 
introduce a simple degenerate spin model and the as- 
sociated dynamics within the degenerate ground states. 
Section IIIII gives a heuristic motivation based on a mi- 
croscopic model more appropriate to real kagome com- 
pounds. In section IIVI we present the results of Monte 



Carlo simulations of the dynamics of the degenerate 
model. In section |Vl we study how the phase space gets 
fragmented in many metastable states and compute the 
configurational entropy from finite-size scaling. We com- 
pare with experiments on kagome compounds in section 
IVII and conclude in section FVlIl 



II. MODEL 

We consider a classical three-coloring model^ with 
spin variables 5*^ = A, B, C (three possible colors, or 
spins at 120°) defined on a lattice, i are the bonds of 
the two-dimensional hexagonal lattice, or the sites of the 
kagome lattice (Fig.[T]). There is a strict local constraint 
which forces neighboring sites to be in different colors, 
and each state p that satisfies the constraint has energy. 
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by definition. The number of degenerate states is macro- 
scopic (extensive entropy) and was calculated exactly in 
the thermodynamic limiti^ As a consequence of Eq. [l] 
the temperature has no effect on the thermodynamics of 
the model: at equilibrium, each state p has the same 
probability. Yet the spin-spin correlations averaged over 
the uniform ensemble arc nontrivial because of the lo- 
cal constraint and decay algebraically ( "critical" state) i^ 
However, the model has no dynamics and one has to spec- 
ify a particular model to study dynamical properties. 

Here we consider the simplest dynamics within the 
degenerate states, i.e. compatible with the constraint. 
While the constraint forbids single color changes, the 
simplest motion consists of exchanging two colors along 
a closed loop of L sites (Fig.[T]). We assume an activation 
process over a barrier of energy kL (where k depends on 
microscopic details), with a time-scale. 



TLiT) = To cxp (kL/T) 



(2) 



where T is the temperature and tq a microscopic time. 
The exact form (Eq. ^ is unessential, the important 
point being that longer loops take longer time (local dy- 
namics). Since the system is known to have a power-law 
distribution of loop length s ^^i^^'^^ (refiecting the critical- 
ity of the thermodynamical state), we have therefore a 
broad distribution of time-scales in the problem. How- 
ever, the loops are strongly correlated and the spin dy- 
namics is nontrivial. 

It has been argued that such constrained problems can 
be described at large scales by effective gauge theories. 
Such examples arc spin-ice systems or hard-core dimers 
which can be viewed as artificial Coulomb phases J^i^ 
The local constraint is solved by an auxiliary (divergence- 
free) gauge field and a long-wavelength free-energy is 
postulated. It describes, as in standard electrostatics, 
algebraic correlations at long distance. The hydrody- 
namic parameters are then extracted from the compari- 
son with exact results (in the present case?^ the Baxter 




FIG. 1: Simplest motion compatible with the constraint: 
color exchange along a loop of length L {e.g. L — 6). Note 
that the flip of the central loop (on the left: before the move) 
facilitates the motion of neighboring sites by creating a new 
fiippable loop (on the right, after the move). 



solution^) or numerics. Furthermore, it also allows one 
to predict a relaxational dynamics {e.g. Langevin) and 
the slowest spin-spin correlations are expected to decay 
as a power-law, as in dimer modelsi^ 

However, we also find a different "short-time" regime, 
resulting from the microscopic model we are considering. 
Indeed, the motion of a loop reorganizes its immediate 
vicinity and can facilitate the motion of a so far frozen 
neighbor (see Fig.[T]). In this sense, this resembles kincti- 
cally constrained models where the motion of a local vari- 
able needs a specific configuration of its neighbors j^^i^ 
but the kinetic constraints here result directly from the 
local correlations. Although the system is fully packed 
with loops (each site belongs to two loops), the issue 
is how the loops (and especially the small loops) self- 
organize. 



III. MICROSCOPIC ORIGIN OF THE MODEL 

We give some heuristic justifications for the model of 
section |lTl based on microscopic considerations. The 
model can indeed be viewed as an effective model within 
the ground state manifold of some more general Hamilto- 
nian, at T ^ JS^ where J is defined below. We consider 
first a Hciscnbcrg model. 



H 



J 2_^ Si • Sj 



(3) 
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where Si is a quantum spin S operator on site i of the 
kagome lattice, and J an antifcrromagnctic coupling be- 
tween nearest neighbor spins. We will discuss the semi- 
classical treatment for which the classical states are the 
important starting point. 



A. Degenerate three-color states 

The minimization of the classical energy associated 
with [Eq. [3] leads to many degenerate states where spins 



point at 120° apart on each triangle. These states are 
not necessarily coplanar; however, the coplanar states 
have the lowest free-energy at low T, a form of (par- 
tial) order-by-disorder to a "nematic" statei^"— Simi- 
larly, for quantum fluctuations at order l/S* in spin- wave 
theory, the zero-point energy is minimized by the copla- 
nar statcsjiii^ However, all coplanar states remain de- 
generate at the harmonic level. The spins pointing at 
120° in the common plane are represented by three col- 
ors A, B, C and the three-color states therefore form the 
ground state manifold of the model [Eq. [3] . 

It is a rather difficult issue to calculate the lifting of 
the degeneracy due to anharmonic fluctuations. In this 
respect, the long-range ordered Neel state with a \/3x \/3 
unit-cell plays a special role. It was indeed argued that 
small-amplitude fluctuations (albeit anharmonic, i.e. at 
the next order in spin- wave theory) favor this state^^"— 
This is similar to the result of Schwingcr-boson mean- 
field theory;^ although this is true only at (small) finite 
TjM From high-T series expansion, the degeneracy is in- 
deed lifted but is a small effect 4i 

By Eq. [I] we assume that the lifting of the degeneracy 
is small compared with both the temperature and the 
energy barriers. 



B. Activation energy, quantum tunneling 

The generation of an energy barrier by fluctuations 
is typical of order-by-disorderj ^^i^^ A canonical exam- 
ple is the Ji — J2 model on the square lattice. While 
two sublattice Neel order-parameters can point in any 
direction at the classical level, the fluctuations select 
the collinear arrangements 1^ The rotation of one sub- 
lattice order-parameter with respect to the other costs a 
(fluctuation-induced) macroscopic energy. There remains 
only two degenerate states separated by an energy bar- 
rier of 0{N), the number of sites (broken symmetry). In 
systems with a macroscopic number of degenerate states, 
the situation is different because local modes connect dif- 
ferent degenerate states. The states are separated by 
barriers of 0(1) and the associated dynamics which con- 
sists of large-amplitude motion of collective spins Eq. [2] 
may be relevant 4^i^i^ There are two different processes: 
the small fiuctuations about a given state of the manifold, 
and the large-amplitude motion within the manifold. For 
continuous spins, the large-amplitude motion consists of 
rotating collectively the spins of a loop, out-of-plane, by 
an angle in a cone at 120°, thus preserving the con- 
straint. The corresponding fiuctuation- induced barriers 
were calculated numerically and appear not to be a pure 
function of the loop length as assumed in Eq. [5J but 
also depend on the configuration;^ However, for small 
loops at the lowest T (when the fluctuation energy is 
dominated by the quantum zero-point motion), E ~ kL 
{k = 0. 14 J5')^ and Eq. [2]is justiflcd. 

At very low temperatures, quantum tunneling through 
the barrier may take placed and the time-scales of Eq. [5] 



saturatei^ The time-scales then depend on the barrier 
shapes and the model consideredi^ 

In real systems, symmetry-breaking fields of spin-orbit 
origin may be present and provide also some energy bar- 
riers. Consider for example, 



H' = H + Dy(S^ 



Z^ 



Ek{<ik ■ Si 



(4) 



which is chosen to be compatible with the 3-coloring 
states: D > is an easy-plane (xy) anisotropy and 
the three vectors dfc are directed at 120° in the kagome 
planci^ In the limit of strong D, H' is analogous to the 
6-state clock model^i except for the degeneracy of the 
classical ground states. We note that in the opposite 
limit of Ising-like XXZ anisotropy, although the system 
orders ferromagnetically, there is a slow persistent dy- 
namics of creation of loops»22 We will restrict the discus- 
sion to £> > and Ek =0 vn the following. 

When the anisotropy is small (which is generally the 
case of intermetallic magnetic ions), rotating the spins 
of a loop continuously by 0, defines a classical energy 
barrier, KLsin^ 9 with k ~ 3DS^/2. 

When the anisotropy is strong (possibly more appro- 
priate to rare-earth compounds), it is too costly to ro- 
tate all components out-of-plane. The lowest-energy ex- 
citations consist of violating the constraint by nucleating 
defects. The simplest effective process for the spins to 
move in the constrained manifold is to create two de- 
fects along a loop. This costs twice the exchange en- 
ergy but the defects are then free to move along the loop 
(deconfinement) and leave behind them a string of ex- 
changed colors.— When the defects annihilate, the loop 
has flipped. The time-scale of this process is given by 
Eq. [21 with K ^ JS^^ This is the important effective 
process in spin-ice systems in general, and the nonequi- 
librium dynamics of defects has been directly studied 
recently<Si24 

Note that by using the discrete model, we intend to 
describe only the slow collective degrees of freedom. The 
rapid motion about the "equilibrium" state (spin-waves) 
is present in the continuous spin model but is integrated 
out in the discrete model (in the barrier), at the first 
order of spin- wave theoryi^ 



to probe the equilibrium state (which is known) but to 
study how the spin dynamics slows down when longer 
loops have to pass higher energy barriers, which take 
more time. The issue is rather to study the relaxation to 
equilibrium. 

The algorithm is similar to that used earlier: (a) we 
choose a single site at random, (b) we choose a neigh- 
bor of this site at random (this defines two colors, hence 
one of the three types of loop A-B, A-C, or B-C) (c) we 
search among its four neighbors the site with the same 
color as the original site (but distinct from it) and we 
iterate until a closed loop is formed (this is guaranteed 
by the periodic boundary conditions) . Contrary to previ- 
ous studies, however, the colors are exchanged along the 
two-color loop (the loop is "flipped") according to the 
probability to cross the barrier, l/r/,. This amounts to 
choose in the Metropolis acceptation rate a microscopic 
rate which depends on the degree of freedom that moves. 
The cluster sizes are N = 3L^, L is the linear size (up to 
L = 144) and a Monte Carlo sweep (MCS) corresponds 
to N attempted updates. 

We have computed the autocorrelation function 



/ 1 "" 



(5) 



where (• • • ) is an average over initial states (10'^ in Fig.[2J 
and up to 10** for better statistics). Here, because we 
have a three-color model, Si(i).Si(0) = 1 for parallel 
spins (same color) and —1/2 for spins at 120° (differ- 
ent colors). By definition, C(0) = 1 and if the state at 
time t is decorrelated from that at t = 0, each spin is 
in one of the three possible colors with probability 1/3 
and C{t) ~ {C{t) measures how long the system keeps 
memory of its initial state). To accelerate the simula- 
tions, we rescale Eq. [2]by Tp = T(,{T) so that the shortest 
loops (hexagons) flip at each attempt. In the following, 
the MCS are in units of t^ and T in units of k. The 
Fourier transform of C(i), C(a;), is the local spin suscep- 
tibility, as measured by experimental probes, for instance 
neutron inelastic scattering (cross-section integrated over 
all wave- vectors), NMR or /itSR on different time-scales. 



A. Summary of the results 



IV. STOCHASTIC SPIN DYNAMICS 

The system evolves in the classical degenerate mani- 
fold by the motion of closed loops, described by a local 
stochastic activated process given by Eq. [2] 

We have studied the spin (color) dynamics by classi- 
cal Monte Carlo simulations. Such Monte Carlo simu- 
lations have been used to study the equilibrium state of 
constrained or loop models,—"— and also in the present 
context.— 1^1^ In these simulations, the updates were 
accepted following the Metropolis algorithm, and irre- 
spective of the length of the loop. Here, the aim is not 



The autocorrelation is given in Fig.[2]for different tem- 
peratures. The relaxation of the system occurs on a time- 
scale Ta and follows a power-law decay, t~'^/'^ (inset of 
Fig. [5]), which is well described by a long wavelength field 
theory, as we shall see. 

Below a crossover temperature T^, the spin dynamics 
develops two distinct time-scales, Tq and Tp {tq, and Tp 
are the notations in supercooled liquids for the long and 
short relaxation times): the autocorrelation decreases 
first into a plateau (quasi-stationary state) and then re- 
laxes to equilibrium. At short times ^ r^, the relax- 
ation is approximately a stretched exponential C{t) w 



exp(— t^) (/? « 0.63). While the dynamics is spatiahy 
homogeneous above T^, it becomes heterogeneous below 
Td with slow and fast regions. 
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FIG. 2: Spin autocorrelation as a function of time (Monte 
Carlo sweeps) with decreasing T (from left to right), C(0) — 1 
{L — 144). Inset: long-time tail (rescaled), 1/f^" (solid line), 
as described by a height model. 



shortest loops. Such definitions and spontaneous genera- 
tions of two time-scales appeared in a different spin model 
where the frustration is played by long-range interactions 
which fragment the system into domains j^ 

By rcscaling all the curves by t^ , we find that the decay 
at long times is a power-law, 



C{t) - 1/ti 



(6) 



with a « 0.33 (see the inset of Fig. [2]). Since a > 0, the 
integrated relaxation time L C{i)dt diverges, and, at 
small frequencies, the Fourier transform diverges like w"" 
(we do not discuss here some natural cutoffs provided by 
e.g. defects at finite temperatures). 

The long-time regime reflects the criticality of the equi- 
librium state and is well described by a free vector-field 
model. The model is obtained by a mapping of the color 
variables onto an auxiliary two-component height field tp 
defined at the centers of the hexagonSi ^^i^^'^" The con- 
struction is as follows: the height vector Lp picks up a 
Ci vector each time it crosses a i =A,B,C color with the 
condition e^i -(- e^ + ec = 0. In such a way, the local 
constraint is automatically satisfied. One assumes that 
the free-energy (of purely entropic origin) reads 



1 



F/T - -K / d^x(V^) 



(7) 
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G-OC(T„)=0.1 
-- 0.42exp(4/T) 
G-O C(T)=0.6 




where is the coarse-grained height field. The stiff- 
ness K = 27r/3 is chosen such as to reproduce the exact 
critical exponent of the spin-spin algebraic correlations, 
■q = 4/3 .^^1^^'^° Eq. ([7]) describes a classical^ interface in 
two spatial dimensions. Similarly to dimer models,— the 
classical fiuctuations of the interface can be described by 
Langevin equations. 
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rf(x,i) 



(8) 



FIG. 3: Two relaxation time-scales. 



where ^(x, t) is a two-dimensional white noise, 
(77(x, i)-'T(x', t')) = TS{x - x')S{t - t'). Eq. [S] describes 
a simple diffusion of the height of the interface. The 
mapping to the slowest spin fluctuations, TOs(x, i) = 
gjQ.<^(x,t)^ IQI = 47r/V3^i2^i^ gives the spin correlations 
at long times and long distance. 



B. Long-time relaxation 

We define the relaxation time of the system, Tq, by 
e.g. C{Ta) ~ 0.1 (the value chosen has no consequence 
as long as it is small enough). In Fig.[3l we give Ta/rp as a 
function of temperature. This ratio becomes much larger 
than one in the limit of low-T, TaJT^ k, 0.42exp(4/T) so 
that Tq ~ '''10 C?") is controlled by the second shortest 
loops (of length 10). For comparison, the time that char- 
acterizes the initial decay of C(i), defined by C(t) = 0.6, 
is of order r^ = rg(r) (Fig. |3]), i.e. controlled by the 



C(x,t) = (m,(x,t)m,(0,0)) 



l-a •' I lljz 



(9) 



with 1 - a = ry/z, z = 2 (from Eq. H]) and /(O) = 1. We 
therefore obtain a = 1/3, in good agreeinent with the 
l/t^~" _ l/f2/3 found numerically (see inset of Fig. [2|). 
The approach also explains that the exponent docs not 
vary with T because the underlying critical phase is in- 
dependent of T, by definition. 

Eq. |9] characterizes the spin fiuctuations at long times 
(by definition of the coarse-grained free-energy) . At short 
times, however, corrections to Eq. [S] are important and 
lead to a different dynamics, as we now show. 



C. Short-time and plateau below Td 



Below a crossover temperature T^ « 1, a shoulder de- 
velops in C{t) and the relaxation time Ta starts to differ 
from Tp which characterizes the initial decay. C(i) devel- 
ops a plateau which becomes more and more stable when 
T is further lowered. The limiting value of the plateau is 
(see Fig. [2]), 



1 ^ 



0.31. 



(10) 



It gives the averaged frozen moment on time-scales 
shorter than Ta, which we note (Si) sa 0.56. On these 
time-scales, only the hexagons have dynamics: all other 
loops are blocked until Tq ~ tio{T) at which a loop of 
length 10 may flip, and the system leaves the plateau and 
returns to equilibrium. 

When the relaxation time of the system becomes longer 
than the experimental time, t^ « r^xp, the system is out- 
of-equilibrium. This occurs at the glassy-like crossover 
temperature, Tg < T^ (which depends on the typical 
time-scale of the experiment) . From the estimation of Tq , 
we have Tg = 10/ln(re:rp/(0.42ro)) w 0.3 for r^^p = lO^s 
and To ~ 10~^^s. For T < Tg the system is trapped into 
the plateau. Once all fast processes have occurred {i.e. 
after r^) the system is in a quasi-stationary state with 
frozen moment squared q (we reserve the term "Edwards- 
Anderson order-parameter" to true equilibrium phase 
transition). 

We can furthermore calculate g as a function of T. It 
is related to the susceptibility by 



(Sf) - (S.) 
T 



T 



(11) 



The frozen fraction depends logarithmically on time be- 
low Tg (see Fig. m, so that x h^^s a cusp at Tg between a 
high-T paramagnetic susceptibility x = 1/2^ 3'iid a low-T 
time-dependent susceptibility. 

The existence of a frozen moment on average is both 
the consequence of frozen regions (which are purely 
static) and dynamical regions with a finite moment on 
average (because of a recurrent behavior). In Fig. [SJ 
we show the autocorrelation Ci{t) = Si(<).Si(0) on each 
site, at intermediate times in the quasi-stationary state 
(—1/2 is white, 1 is black if it has never moved between 
and t or gray otherwise). While most sites have dy- 
namics (white and gray), there is a fraction of frozen 
sites (in black). The averaged fraction of frozen sites is 
Nf = 0.121A^, and the probability distribution function 
is found to be gaussian (as a consequence. Fig. [5] is typ- 
ical of what happens at low T). The existence of 12.1% 
of frozen sites explains only part of the averaged frozen 
moment, q = 31%. In addition, other (dynamical) sites 
contribute. This is because the frozen clusters provide 
boundary conditions for the neighboring sites and the 
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FIG. 4: Frozen moment squared as a function of temperature 
and observation time, Eq. 1101 L — 144. 



constraint propagates between clusters. For instance, the 
spins on the outer side of the cluster boundary can take 
only two of the three possible states, the third possibility 
being frozen inside the cluster. They have hence stronger 
probabilities to return to the original value. In Fig. [51 we 
indeed see, extending between the frozen clusters, large 
dynamical regions where the spins are in their original 
state (gray). These constrained regions contribute to al- 
most two thirds of the averaged frozen moment. 




FIG. 5: Real space picture of the autocorrelation, Ci{t) = 
Si(i).Si(0), at time t = 10^ and T = 0.1. Black: frozen 
sites. White {Ci{t) = —1/2), in gray the sites which have 
moved between and t but have returned to their initial state 
(a(t) = 1). 



Furthermore, it is seen in Fig. [S] that frozen sites form 
clusters randomly distributed over the system. The num- 
ber of spins in a cluster is distributed according to Fig. [51 
The average is (s) = 42 sites (and is size independent for 
L > 72, sec inset of Fig. 15]), thus defining an emergent 




FIG. 6: Distribution of tlie sizes of the frozen clusters. The 
average is (s) — 42 sites (inset: finite-size effect) or the emer- 
gent lengtli scale is (s)^ . The dashed line is a guide to the 
eye (exponential). 
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FIG. 7: The radial distribution function of active degrees of 
freedom: probability to have a flippable hexagon at distance r 
from a given flippable hexagon at (normalized by the num- 
ber of hexagonal sites) averaged over the uniform ensemble. 
While the nearest neighbor position is not compatible with 
the constraint, the first peak corresponds to an attraction of 
next nearest neighbor hexagons. 



length scale (s)^'"^ = 6.5 intersite spacings. The pic- 
ture of the frozen phase is that of "jammed" clusters of 
nanoscopic scale (s)^^^ occupying 12.1% of the sites. 

What is the origin of the "jamming"? First, "jammed" 
clusters do not contain flippable hexagons (by definition) 
but are, of course, criss-crossed by longer loops which 
are blocked at the temperatures considered. This im- 
plies that a typical three-coloring state must have a low- 
enough density of flippable hexagons. On the kagome 
lattice the density of flippable hexagons (averaged over 



the uniform ensemble) is 0.22, so that forming a large 
cluster of (s) =42 sites on average is unlikely in absence 
of correlations. In Fig. [3 we give the correlations g{r) 
(radial distribution function) in the positions of the flip- 
pable hexagons of the same type<^ We find indeed a 
strong attraction: the neighboring hexagons cannot be 
occupied by the same type of loop (it is incompatible 
with the constraint) but the second neighbor positions 
are highly favored (attraction). There is a high probabil- 
ity to have a fiippable hexagon if the (second) neighbor is 
a flippable hexagon. This attraction creates aggregates 
and voids, opening the way to regions free of flippable 
hexagons. The system can therefore be viewed as a mi- 
croscopic phase separation of active and inactive regions, 
the active regions having flippable hexagons, the inactive 
regions having longer loops. Recall that the degenerate 
model can be seen as being at the boundary of a phase 
transition in parameter space^^ in particular between 
active and inactive phases, having respectively short and 
long loopsi^ This is a necessary but not a sufficient con- 
dition for the region to be "jammed" because the number 
of flippable hexagons is not conserved by the dynamics 
and they "move" on the lattice (see Fig. [T|). The frozen 
clusters correspond to special configurations and regions 
inaccessible to flippable hexagons. For example, a frozen 
cluster of 12 sites is shown in Fig. [5] each hexagon on the 
border has the three possible colors A, B, C, thus making 
it impossible to create a flippable configuration. One can 
have clusters of arbitrary size (see Fig. [5]) or walls that 
prevent flippable hexagons to diffuse in different regions 
of the sample. 




FIG. 8: "Jammed" cluster (the smallest one, in gray): no 6- 
loop can unjam any of its 12 sites. The shortest "unjamming" 
loop (of length 10) is shown (dashed line). 

However, a loop of length 10 (shown by a dashed line 
in Fig. [5]) will unjam the conflguration, and the clus- 
ter shown will be annihilated. The way the relaxation 
takes place at longer times is via the dynamics of cre- 
ation and annihilation of "frozen" clusters on time-scale 
Tio(T). For T < Td, there is a separation of time-scales 
between the "rapid" hexagon motion ~ tq{T) and the 
longer creation/annihilation of frozen clusters ~ rio(T). 



D. Dynamical heterogeneities T < Td 

We now consider some dynamical local quantities. Fol- 
lowing studies of standard glasses;^ we define a local mo- 



bility field Ki{Q,t) which measures how many times the 
site i has changed color during the time interval between 
and t. It is linear in t for large t so that one can define a 
local frequency fi = Ki{0,t)/t. Frozen sites have fi = 
while dynamical sites have ft > 0- 

FIG. 9: Dynamical heterogeneities in space. The gray scale 
is proportional to the local frequency fi of the site from black 
(frozen) to white (high frequency). 
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FIG. 10: Non-uniform slowing down of the dynamics by lower- 
ing the temperature. From T > Td, homogeneous (gaussian) 
distribution of local frequencies to T < T^, heterogeneous 
(skewed) distribution and at T < Tg a frozen fraction ap- 
pears, t = 10* and L = 144. 

The real space picture of /; at a given time is given in 
Fig. m from black (frozen sites) to white (fast sites): we 
see the variations of the local dynamics across the sys- 
tem and some clusters of slow frequencies, i.e. a form of 
dynamical heterogeneity. We plot the corresponding his- 
togram of frequencies in Fig. [TO] at various temperatures. 
At high temperature, the distribution is homogeneous 
(gaussian). At lower temperatures, the dynamics slows 
down, and the distribution broadens and becomes asym- 
metric (non zero third moment or skewncss). Eventually 



at T < Tg, a frozen fraction appears and the distribution 
becomes continuous between two typical peaks. 



^(/) = ^5(/)+^(/) 



(12) 



where A(/) is a smooth broad function. One can de- 
scribe this evolution as a crossover between a homoge- 
neous high-temperature phase with a single type of dy- 
namical site and a low-temperature phase with many in- 
equivalent dynamical sites. It can be described in terms 
of large-deviation functions and a "free-energy" can be 
defined,^ 



V. FRAGMENTATION OF THE PHASE-SPACE 

We show that the phase space is fragmented into a 
gWSc number of sectors for T < Tg, separated by barriers 
of 0(1). For this, we directly enumerate all the states of 
small clusters and analyze how the system evolves in the 
phase-space as a function of temperature. This allows us 
to describe the landscape of energy barriers separating 
states and basins, i.e. a hierarchical organization of the 
states (nonfractal here). 

Let Pp{t) be the probability of the system to be 
in a configuration p ~ 1, . . . , Nc, Nc the total num- 
ber of states which we have numerically enumerated on 
small clusters with periodic boundary conditions (N ~ 
27,36,81,108). We have found Nc = 6.4 x 1.122^^ 
(dashed line in Fig. [TT|) . slightly smaller than the exact 
result in the thermodynamic limit 1.135^i^ 

The master equation governing the dynamical evolu- 
tionofP(t) = (Pi(t),...,PjVcW), 



dP 



(13) 



involves a matrix w which contains the transition rates 
from a configuration ptop'. The only allowed transitions 
are single flips of loops of length L, Wp^p' = —1/tl{T) 
where tl (T) is given by Eq. [2j Here from detailed bal- 
ance, we have Wp^p' = Wp>^p {Ep = for all states) 



and w. 



p-).p 



^ 'p'- 



^p Wp^p' 



ensures the conservation of 



the probability, '^pPpif) = 1. All states satisfying 







(14) 



are stationary, such as, in particular, the equilibrium uni- 
form distribution Pp{t) — 1/Nc. w may have more than 
one zero eigenvalue and the additional stationary states 
prevent the system from exploring the phase space (bro- 
ken crgodicity). Examples are systems with a broken- 
symmetry, the phase space of which has a finite number 
of disconnected sectors in the thermodynamic limit. In 
each sector the Gibbs distribution is stationary assuring 
as many zero eigenvalues as the number of sectors or bro- 
ken symmetries. By contrast, in a glassy-like phase, the 
number of eigenvalues satisfying e <C ^/Texp (jexp is the 



experimental observation time) scales like e 



NSa 



there 



is a finite configurational entropy Sc- In other words, a 
macroscopic number of states, thus differing at the mi- 
croscopic scale, never relax on the observation time-scale. 
In the present model, w has a finite hierarchical struc- 
ture. Here it is a consequence of a microscopic model and 
is not assumed from the beginning as in hierarchically 
constrained modelsj^i^i^ Contrary to these examples 
(or spin glasses^), however, we find only four levels of 
hierarchy: the phase space is split into a few "Kempe" 
classesr^^ which are split into ^ N topological sectors 
and then in e^'^'= trapping sectors (sec Fig.[Tl]for a graph- 
ical illustration of this hierarchy in the phase space). 
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The dynamics of loops of all sizes is known to be noner- 
godic on the kagome latticci^— It means that moving all 
loops is not sufficient to go from a given state to any other 
state in the phase space, w split in "Kempe" classesj^i^ 
the number of which is in general unknown.— 

Since it is therefore impossible to enumerate all states 
by moving loops iteratively, we have allowed to intro- 
duce defects that violate the 3-colored constraint. To 
control the density of defects, we have introduced an en- 
ergy penalty, i.e. the antiferromagnetic three-state Potts 
model. By cooling the system at low temperatures in 
a Monte Carlo simulation, one generates three-coloring 
ground states that are in different "Kempe" sectors (and 
the sectors themselves by switching on the loop dynam- 
ics). For N = 108, we find four sectors, a large one with 
89% of all states and three smaller ones, all separated by 
infinite barriers for the loop model. 

Within each Kempe sector, the three-coloring states 
can be characterized by topological numbers. They are 
defined by counting the number of colors along nonlocal 
horizontal and vertical cutsi^ There are six such num- 
bers w^'^ {i = 1, 2, 3), which may take any integer value 

from to L with the constraint X]i=i '^V^ ~ -^' ^^ ^^^^ 
of them are independent. This gives at most N'^ sectors, 
but since some combinations are not allowed, the num- 
ber is of order A^ (Fig. [TT|) . The dynamics of local loops 
conserves these numbers so that each Kempe sector is di- 
vided into N topological sectors. Only winding loops of 
length L or L^ (the longest loop takes all two color sites 
and has length 2N/3) may change them. In fact, the av- 
eraged length of the winding loops scales like J^^/^ , 48,52 
The topological sectors are therefore separated by barri- 
ers growing with the system size like L^'"^, defining in- 
finite barriers in the thermodynamic limit and broken 
ergodicity sectors. This is analogous to the "jamming" 
transition induced by additional forces: the favored or- 
dered state needs rearrangements of infinite loops in or- 
der to equilibrate.— 1^ Here we recall that the phase space 
is in general broken into ^ N sectors (which we have ex- 
plicitly constructed), labelled by quantities conserved by 
the local dynamics 4^ 



FIG. 11: Hierarchical structure of the phase space and ex- 
ponential number of disconnected sectors Nt at low tem- 
peratures (solid circles). Nt corresponds also the number 
of zero eigenvalues of the matrix w. The total number of 
states Nc (squares), topological sectors Ntopo (diamond) re- 
sult from the complete enumeration of states on clusters of 
size Af = 27,36,81,108. 



B. Fragmentation in e^ "^ sectors 

For T < Tg, the dynamical matrix w split further into 
new smaller sectors which we have constructed for dif- 
ferent system sizes. We find that the phase space is 
split into 1.085^ independent trapping sectors fFig. lTl]) . 
The spin dynamics has a fast equilibration within a sec- 
tor characterized by the motion of 6-loops on time-scale 
Tp = TeiT) and the motion between sectors occurs on 
time-scale, t^ ~ tio(T), which is frozen below Tg by 
definition. Above Tg, the system equilibrates within a 
topological sector. 

The number of sectors defines a finite averaged con- 
figurational entropy per site Sc = In 1.085 ~ 0.082, 
which is approximately two thirds of the full entropy 
Seq = In 1.122 = 0.115. Upon reducing the tempera- 
ture, the system goes from an equilibrated state with the 
full entropy Seq (the number of topological sectors is sub- 
extensive) to a metastable state below Tg where it looses 
the configurational entropy: 



AS = Sc = 0.082 = O.lSeq 



(15) 



The configurational entropy refiects in phase space the 
entropy of the microscopic arrangements of the frozen 
clusters (section ITV) . A crude comparison consists of dis- 
tributing Nf/{s) disks on the lattice {Nf/{s) is the num- 
ber of frozen clusters of average size (s) = 42, we note x 
the density), with entropy S/N ~ (— xlnx— (l — x) ln(l — 
x))/{s) = 0.009 {x = 12%). This is too small though by 
an order of magnitude compared with Sc- 

For Tg < T < Td, one can define coarse-grained 
states by eliminating the fast dynamics into an entropy. 
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While, on average, each sector contains (1.122/1.085)^ — 
1.034^ states (thus defining the averaged entropy 52 = 
0.0347V), we find a broad distribution of sector sizes from 
s = 1 (a single state) to a large sector s < Nc- How- 
ever, we believe that this is a finite-size effect. Indeed the 
probability to fall into a sector of size s is found to be 
roughly constant at small s and increases for larger sec- 
tors. In contrast, for a Monte Carlo sampling of states 
as done in section ITVl the frozen fraction distribution is 
homogeneous (gaussian) for L > 18, while for L < 18, a 
large portion of states has no frozen fraction at all. As 
a consequence, the distribution of entropies is certainly 
more homogeneous for large system size. 

In summary, we find that the phase space has hierar- 
chical levels: it has sectors characterized by conserved 
quantities and separated by infinite barriers (broken er- 
godicity) and sectors or traps separated by finite barriers. 
The number of topologicals sectors is of order N (noncx- 
tensive entropy), and there is no essential difference be- 
tween them at the microscopic or mesoscopic scale: a 
local measurement cannot distinguish between two dif- 
ferent sectors. On the other hand, the number of traps 
is of order e''^'^'^ (finite configurational entropy). There- 
fore the system looses a finite entropy at Tg and a local 
disorder is self-induced: a local measurement can dis- 
tinguish between two mctastable states (for instance, if 
there is or not a frozen cluster). In this sense, Tg can 
be called a glassy crossover temperature. By opposition, 
the jamming transition found in Refs. I48ll49l corresponds 
to broken ergodicity associated with a sub-extensive en- 
tropy (no self-induced disorder). 



VI. DISCUSSION OF EXPERIMENTS 

We now discuss the kagome compounds that have a 
freezing transition. We argue that the freezing temper- 
ature Tg is governed by the energy scale of the barriers 
and when possible, we identify the possible mechanisms 
we have discussed in section IIIII the barriers are either 
dynamically generated by the rapid spin- wave motion or 
generated by anisotropics, depending on specific materi- 
als. We also compare the strength of the "frozen" mo- 
ment to the experiments available and the dynamics of 
the system. Note that the present dynamics of loops 
is classical (if a quantum coherence is maintained, the 
system was predicted to order—). Some quantum fluctu- 
ations are therefore neglected here, but may turn out to 
be important, especially for the copper oxides discussed 
below {S = 1/2), if the anisotropy is small enough4^ 



A. SrCrgpGaia-gpOig (SCGO) 



In SCGO, a phase transition occurs at Tg ~ 3.5 — 7K, 
depending weakly on the Cr'^+ {S = 3/2) coverage pJ^"— 
Tg depends also on the experiment: Tg ~ 3.5K by squid, 
5.2K by neutron scattering for the same compoundi^ 



What could be the appropriate microscopic model? 
The Cr'^'^ ions have no orbital moment {L = 0) and 
the spin anisotropy is expected to be small. From 
EPR indeed, DS^ ^ 0.2K,^ By constrast, the mea- 
surements of the spin susceptibility on single crystals 
showed a large anisotropy disappearing when increas- 
ing the tcmperaturci^ This was therefore attributed to 
the spontaneous breaking of the rotation symmetry by a 
nematic order (coplanarity), and not a real anisotropy 
of the modeli^ Similarly, the 8K barrier obtained by 
/iSR for p — >■ which was originally interpreted as a 
large single-ion anisotropy?^ is in fact absent if one uses 
a different fit of the datai^ On the other hand, for 
p — ^ 1, energy barriers of '--^30K were obtainedi ^^'^'^ Since 
they are two orders of magnitude larger than the spin 
anisotropy, they are more likely to be induced by the 
fluctuations. With E = kL = 30K and L = 6, we 
have Tg = 0.3k = 1.5K. On the other hand, if we use 
K = 0.1475 (section IIII|) and J ^ 50K from the spin 
susceptibility, we find Tg = 0.0475 ~ 3K. Both esti- 
mates are in fair agreement with the experimental result. 
The model does not predict, however, a thermodynamic 
transition while, experimentally, this has been a disputed 
point, especially regarding the sharpness of the nonlin- 
ear susceptibility xs'^'^ We also note that the "thermo- 
dynamic" anomalies we have mentionned at Tg are not 
only rounded but also the entropy change AS = 0.082Af 
is small compared with the full entropy A^ ln(25 + 1) of 
continuous spins. Yet this makes a definite prediction for 
the entropy change. 

Furthermore, the frozen moment measured in neutron 
elastic scattering is small, (5;)^ ~ 0.12 — 0.24 of the max- 
imum moment (depending on the Cr coverage) and most 
of the signal is in the inelastic channeh ^^i^^'^" In the ex- 
perimental setup of Ref. [13, the inelastic channel starts 
above the neutron energy resolution of 0.2 meV, giving 
in that case a lifetime of the frozen moment longer than 
~ 20ps. Neutron spin-echo showed that the moment is 
still frozen on the nanosecond time-scale at 1.5 K3i How- 
ever, no static moment was originally observed in ^SR,— 
but a weak static component may not be excludedi^ Sim- 
ilarly, in Ga NMR. the wipeout of the signal shows a 
dynamics that has slowed down but is still persistent i^ 
However, in both cases the muon or the Ga nuclei probe 
many sites and may see primarily the dynamical sites. 

In the model developed above, the system remains dy- 
namical below Tg. The system has flippable hexagons 
on time-scale rg (T) but also spin- waves on a more rapid 
time-scale, which we have not described. The latter 
should contribute to the specific heat as in normal two- 
dimensional antiferromagnets and give in particular a 
T^ specific heat as observed experimentallyfi^ This is 
the consequence of the two Goldstone modes associated 
with the selection of a common plane (nematic broken 
symmetry) .— 

We can make different assumptions regarding the time- 
scale of the activated dynamics with respect to the ob- 
servation time-scale. If tq{T) ^ Tneut., the system is 
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trapped into a typical 3-coloring on the experimental 
time-scale. Still the averaged moment is different from S 
because of the rapid zero-point fluctuations of the spin- 
waves. One can estimate that the effect of the two Gold- 
stone modes is to reduce the moment to m = S — 0.16i^ 
For Cr'^^ [S ~ 3/2), the correction is small and cannot 
explain the small moment measured. 

Suppose now that the hexagons still have a dynamics, 
as indeed predicted for T < Tg. We found in this case 
that the frozen moment is (Si)'^ f« 0.31 (Fig. [2]). Apply- 
ing the same zero-point motion reduction as above, we 
find 0.31(1 - 0.16/S)^ = 0.25 which is close to the ex- 
perimental frozen moment. The model therefore gives a 
fair account of the measured frozen moment. The small 
static moment is not due to strong quantum fluctuations 
but rather to the loop (hexagon) fluctuations. 
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FIG. 12: Spectral function at diflterent T (dashed lines are 

To characterize the dynamics, we have computed the 
local dynamical response at different T fFig. [T2|) . These 
are the Fourier transforms of the autocorrelation func- 
tions given in Fig. [5J At T > Tg, and low frequencies, 
we have C(w) ^ w"^/'^ as the consequence of the univer- 
sality of the height model. However, this is valid over a 
limited range of frequencies: in Fig. [T^l the dashed lines 
give examples of power-laws with exponents 0.33 and 0.7, 
for comparison (note that all the curves are shifted hor- 
izontally by 1/t^). It is also in fairly good agreement 
with the observed power law behavior in neutron inelastic 
scattering on powders, w""'* above the transitiouji^r— 
When T is lowered, the quasi-elastic peak corresponding 
to the frozen moment develops. Note that the sum-rule 
J C{uj)duj = 1 ensures that the apparent loss of intensity 
at low temperatures in Fig. [T^] corresponds to a transfer 
into the elastic peak. Although the approach is different, 
we note that the exponent is not far from that obtained 
by dynamical mean- field theory, a ~ 0.5i^ 

In summary, the model describes a dynamical freezing 



crossover into a partially frozen phase and a small frozen 
moment in overall agreement with the experiments. The 
broad neutron response is interpreted as the motion of 
loops above Tg. In the frozen phase, only the hexagons 
are predicted to move (in addition to spin waves). They 
could possibly be characterized by special magnetic form 
factors, as in ZnCr204i^ 



B. Volborthite Cu3V207(OH)2.2H20 

In volborthite j2^ a freezing transition occurs at Tg ^ 
IK, with a finite static moment observed by NM R^^i^^ 
but no long-range correlations in neutron scattering ii2£ 
Volborthite is a slightly distorted kagome lattice and 
there is some current debate as to whether the main 
magnetic couplings arc kagome like or more one- 
dimensional»^ We will assume below that it can be 
viewed as a kagome antiferromagnet and that the dis- 
tortion is a small effect. 

Below the transition, NMR revealed that the phase is 
heterogeneous with a time-dependent lineshapc, leading 
to distinguish between "fast" and "slow" (static) sites, 
either at small fielda^S or in a dinstinct phase?^ at larger 
fields,^ 

These results resemble the dynamical heterogeneities 
found in the model below Tg. We can make a more de- 
tailed comparison by computing the distribution of fields. 
NMR was performed on vanadium nuclei which seat at 
the centers of the hexagons i ^^i^^ The nuclei see effective 
fields averaged over the six sites in of a hexagon H (as- 
suming for simplicity the same hyperfine coupling Ai^ ) , 



{hH)= Y.A.^H- I dt%,it'), 

A ,, 1 "^ 



(16) 



which depend on the hexagon (inhomogeneous broaden- 
ing) . An average over the NMR time scale t is taken. In 
principle t is much larger « 10 — 100/is than the micro- 
scopic time-scales ~ps, and t can be taken to -|-cxd. In 
systems with slow dynamics, NMR probes local trajec- 
tories averaged over t. The lineshape depends on i, thus 
providing information on the presence of dynamical het- 
erogeneities. The lineshape is related to the distribution 
function of field strengths P{h = |(h//)|), which we have 
calculated in the present case. 

We expect different regimes, according to whether the 
NMR time-scale t is shorter or longer than the charac- 
teristic time-scales of the dynamics, r^ and Tq. Note 
that since these describe activated processes, they may 
become much longer than the ps microscopic time at low 
temperatures. 

• t ^ Ta,Ti3. The system equilibrates on NMR time- 
scales, e.g. at high-T. Every site has dynam- 
ics and summing random vectors (at 120° though) 
gives a gaussian distribution of fields (dashed line 
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FIG. 13: Distribution of local field strengths at the centers of 
the hexagons. Evolution from high T to T < Tg. A typical 
static state (no loop flip) is given for comparison (peaks at 
0, \/3, 3, see top). 



in Fig. [13]). For t —> oo, summing local fields corre- 
sponds to a random walk and the typical strength 
h ^ 1/Vt — > since we have no external field. 

• t <C Ta,Tp. The system is completely frozen in 
a typical three-coloring state. Each nucleus sees 
a well defined static field. For a 3-coloring, there 
are only three possible field strengths at the center 
of the hexagon, /i = 0, VS, 3 (see the configura- 
tions shown in Fig. 1131 top). Averaging over the 
uniform ensemble, we find three peaks with weight 
18%, 60%, 22% (22% is the fraction of flippablc 
hexagons). For comparison, the Q = Q antiferro- 
magnetic state would have a single peak aX h — Q 
with 100% of the hexagons and the V3 x V3 state 
a single peak at /i = 3. 

• T;3 <C i <C Td . The system is out-of-equilibrium be- 
low Tg, by definition. The dynamical sites provide 
a time-dependent averaged field (broad part of the 
lineshape in Fig. [T3|) . The frozen sites inside the 
clusters provide a static field: we find two peaks 
at /i = and v3 and no peak at /i = 3 which 
corresponds to the flippablc hexagons. Although 
the local field does not change when they fiip, the 
probability that they remain in a fiippable config- 
uration is small. Instead they move on the lattice 
and there are very few isolated flippablc hexagons 
inside frozen clusters. We further note that the 
static fields inside the frozen clusters show a ra- 



tio P(0)/P(V3) sa 0.8 much larger than that of a 
typical state « 0.3 (Fig. [T5)) . This means that the 
frozen clusters resemble locally to the (5 = state, 
the state with long linear winding loops, precisely 
those which do not fiip. 

Experimentally, in volborthite, the NMR lineshape 
consists of two dynamically heterogeneous contributions 
at T < Tg 1^ A slow "rectangular" shape was assigned to 
a static field (the "rectangular" shape arising from the 
powder convolution) and a fast gaussian contribution to 
dynamical sitesi^ Similar results were obtained at low 
fields on lower-quality samples i^ 

The low-field results can be compared with Fig. [13] 
(up to a powder convolution). For T < Tg we find two 
contributions: (i) a static contribution coming from the 
frozen clusters and represented by two peaks at h = and 
h = v3, the latter giving rise to a rectangular shape in a 
powder sample, (ii) a dynamical part resulting from the 
dynamical regions and giving rise to a broad response. 
Here we do not have a single dynamical site but this is 
not necessarily incompatible with the experiments be- 
cause of the difficulty of resolving different sites. We also 
note that the respective contribution of both is smaller 
than the 50%-50% observed experimentally;^ but this 
depends on the time-scale. 

Assuming that the system had a -y/S x -\/3 character and 
that the static field was therefore due to the h = 3 types 
of hexagons, a small frozen moment of OAlfis per site 
was extractedi^ In the present model, a peak at ft, = 3 
is not compatible with the existence of dynamical sites. 
Instead we assign the experimental peak to the h = v3 
frozen field. In this case, instead of m = 0.41^b^ the 
static moment is tti = 0.41^b x \/3 = 0.71/x_b (as also 
proposed in Ref. |43 for different reasons), which is more 
compatible with conventional on-site zero-point fluctua- 
tions, (1 — 0.16/S')/i_B — 0.68/i_B. Moreover, if we now 
calculate the total frozen moment averaged over all sites 
as measured by neutrons (while NMR sees the full local 
frozen moment), we would predict rriav = 0.56x0.71/iB = 

0.40^B- 

We conclude that the present study gives a model 
for the phase transition and the heterogeneous state ob- 
served in volborthite. Similarly to SCGO, it gives an in- 
terpretation for the small moment observed for T < Tg-. 
the fluctuations of small loops reduce the averaged mo- 
ment. The model suggests a more precise picture of 
frozen clusters with an emergent length scale, that can 
be further tested experimentally. 



C. Vesignieite Cu3BaV208(OH)2 

For vesignieite^ Tg — 9K, and the ground state is also 
heterogeneous: approximately 50% of the sites (muon 
sites and nucleus sites) experience a static field.— i^ The 
loss of 50% of the total intensity in NMR is due to the 
nuclei which have a time-scale that cannot be detected, 
and therefore reflects some dynamical heterogeneities in 
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the local environment. It would be inaccurate to consider 
that the 50% of the observed intensity is due to the spins 
in the frozen clusters and the missing 50% due to the 
fast moving spins. It may well be that some dynamical 
sites of Fig. [13] are detected (this is in fact what we as- 
sumed for the volborthite where 100% of the nuclei were 
detected). The fact that the fraction docs not match the 
number of frozen sites of 12% is not therefore a serious 
drawback. Alternatively the fraction of frozen sites cer- 
tainly depends on the interactions. For Dzyaloshinskii- 
Moriya interactions which are present and may be rather 
strong;^ the frozen fraction will certainly increase be- 
cause it favors the Q = state with long loops. 



D. Hydronium jarosite, (H30)Fe3(S04)2(OH)6 

A freezing transition occurs at Tg ^15Ki^ By vary- 
ing sample preparations, Tg ~ 12K-18K and it appears 
to be weakly sensitive to the Fe coverage in the range 
92-100%!^ Neutron elastic scattering has found short- 
range correlations of the -^/S x v^ type, but no long-range 
order 4aiiHii2^ The Heisenberg coupling is JS^ == 244K 
(5 ^ 5/2)^ so that Tg/JS^ = 0.05 which in terms of a 
classical Heisenberg model means that the system should 
be in the collective paramagnetic regimei^ One clearly 
needs some additional ingredients to explain the freezing 
transition. 

Spin anisotropy is known to be present in a similar 
jarosite compound, KFe3(S04)2(OH)6, both single-ion 
easy-plane anisotropy DS^ ~ 30K and a Dzyaloshinskii- 
Moriya interaction IZJIS*^ ^ 20K explaining the excita- 
tion spectrumii2iii£^ X-ray dichroism of the Fe'^"'" ion 
also found a single-ion anisotropy in good agreement with 
the above figureJ^ In addition, in ordered jarosite com- 
pounds, a second transition corresponding to the in-plane 
locking of the spins occurs at 45 — 55Kji22 With these 
large values in mind, we assume that the energy barriers 
of the model originate in the anisotropy. In this case, 
we can predict Tg and compare with that obtained from 
ac-susceptibility measurements i^ Since, in the model, we 
have Tg « 0.3k = 0.225DS'^ for r^^p = lO^s (6x lO-^Hz), 
and DS^ ~ 30 - 55K, we find Tg « 7 - 12K. Similarly, 
for Teoop = 80ms (80Hz), we find Tg « 9 - 15K. These 
estimates are a little smaller than the experimental fig- 
ures and depends more strongly on the measurement fre- 
quency (the same distinction occurs in structural glasses 
between "fragile" and "strong" glasses) i^ Moreover, by 
varying synthesis conditions, Tg was found to be cor- 
related with the distortion of the FeOe octahedra: the 
stronger the distortion the larger the Tg <^^ Since the oc- 
tahedron distortion implies a linear change in the crystal 
field splitting, hence in the single-ion anisotropy D, we 
expect indeed linear changes in Tg w D, as observed 
exp er iment allyJ^ 

For T < Tg, an estimate of the frozen moment has 
been obtained by /iSR and amounts to 3.4/iB compared 
with 5.92^B of the Fe3+ ion^ so that (S^) = 0.57. It is 



not far from the present estimate 0.56(1 — 0.16/5) = 0.52. 
However, it is surprising that similar values were obtained 
in ordered jarosites J^ 

For T > Tg, neutron inelastic scattering has been per- 
formed and showed that the local response, x"(a;) 



„ -0-68^ ^-j- very low frequency, we have found uj"^'^ 
but at larger frequencies it could be fitted by a larger 
exponent (the second dashed line in Fig. [T^] corresponds 
to a = 0.7). The agreement is therefore qualitative with 
a broad increasing response by lowering the frequency 
(to be contrasted with the flat response of a conventional 
two-dimensional antifcrromagnct) but a single exponent 
is not found. 

To conclude, the present study suggests that Tg in 
(H30)Fe3(S04)2(OH)e is related to a dynamical freez- 
ing into a heterogeneous state. The relevant energy scale 
here, contrary to SCGO, is the anisotropy, as experi- 
mentally claimed ji2^ Below Tg, we expect a small frozen 
moment on average and a persistent dynamics of the 
hexagons, which distinguishes the present transition from 
a complete dynamical arrest. More studies of the low- 
temperature phase would be interesting. 



E. Other kagome compounds, competitions 



It is well known that not all kagome compounds 
have a freezing transition and we briefly discuss some 
other compounds. Some have magnetic long-range or- 
der, which is often accounted by additional spin interac- 
tions. Others, such as the herbertsmithite compounds 
ZnCu3(OH)6Cl2iii and MgCu3(OH)6Cl2^ have no 
freezing transition (unless an external field is appliedi^) 
and no long-range orderJ^ The neutron inelastic re- 
sponse has no clear energy scale in ZnCu3(OH)6Cl2^ 
and is fitted by a broad power law u!~'^'^'^ at low enough 
energy)^2i^ with some similarity with that of SCGO and 
the hydronium jarosite above Tg. In the present model, 
one would interpret this result as being in the phase 
above Tg , and the neutron inelastic response agrees qual- 
itatively with Fig. 1121 However, the reason why Tg would 
be smaller than the lowest temperatures reached exper- 
imentally, say 50 mK, is not clear. We have argued 
that Tg is controlled by the anisotropy (dynamically- 
generated or not), and the anisotropy is present in 
ZnCu3(OH)gCl2J ^^"'^'^^^ Two important effects are miss- 
ing: it is known that antisite disorder is present?^^ 
and that S = 1/2 compounds have strong quantum ef- 
fects with currently debated quantum spin liquid phases 
if the anisotropy is sufficiently weak (such a coupling 
may discriminate between different phases m S = 1/2 
compounds^). It is therefore clear that competitions 
arc important to account for all these phases. 
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VII. CONCLUSION 

We have described a simple spin model which has a 
dynamical glassy-like freezing at a crossover temperature 
Tg, in absence of any quenched disorder. 

The system evolves from a dynamically homogeneous 
phase with a single time-scale (T > Td) to a dynamically 
heterogeneous phase with two time-scales (T < Td). The 
first time-scale r^ ^ ''6(2^) corresponds to the "rapid" 
degrees of freedom, the shortest loops. The second time- 
scale Ta is associated with the rearrangement of the 
"frozen" clusters. The frozen clusters have a microscopic 
length-scale (they typically contain a few tens of sites) 
but their rearrangement time is not controlled by the 
their size but by the size of the second shortest loops, 
Ta ~ tiq{T). When Tq becomes longer than the ex- 
perimental time-scale for T < Tg, the system is out-of- 
equilibrium and glassy-like. The clusters contain spins 
that are frozen on the experimental time-scale and real- 
ize a microscopic-scale disorder. In this case, the system 
has a finite (small) averaged frozen moment but no true 
long-range order. We have explained that the frozen mo- 
ment is due, partly to the frozen clusters themselves and, 
partly to dynamical regions where the spins are strongly 
constrained by the frozen regions. 

The phase space of the system appears to be orga- 
nized in a partially hierarchical manner with conserved 
quantities defining ~ N basins separated by infinite bar- 
riers (broken ergodicity). Each basin was shown to fur- 
ther split into e'^^" sectors separated by finite barriers 
which trap the system in a metastable state below Tg. 
This macroscopic fragmentation of the phase space cor- 
responds to the local disorder induced by the "frozen" 
clusters. At Tg, the system has therefore some "thermo- 
dynamic" anomalies characterized by the loss of the con- 
figurational entropy which we have calculated by finite- 
size scaling, Sc = 0.082 per site. 

The system undergoes a glassy-like transition at Tg be- 
cause the residual "rapid" degrees of freedom (the short- 
est loops) only partially reorganize the system. In a typ- 
ical state, the density of the shortest loops is not very 
small, but, by effectively attracting each other, they form 
aggregates and voids (micro phase separation), the latter 
regions being, hence, frozen. Some details, as to what 
their density is, or how they precisely interact, certainly 
depend on the system and the model, but the mechanism 
we have presented here is rather clear: the strong local 
correlations generate slow extended degrees of freedom 
which, since they are correlated and attract each other, 
"phase-separate" in dense active regions and void inac- 
tive regions. 

Several aspects of the degenerate model are simply as- 
sumed. We have assumed the absence of long-range order 
by considering degenerate states (Eq. [T]), and an acti- 
vated relaxation time (Eq. [2|); hence, not surprisingly, 
the dynamics is slow. We have discussed in section IIIII 
why both assumptions may be approximately realized in 
microscopic models with continuous degrees of freedom. 



We argued that the origin of the energy barriers is the 
partial order-by-disorder, i.e. the barriers are dynami- 
cally generated by the rapid spin- waves, or by an explicit 
anisotropy arising from the spin-orbit coupling. The de- 
generacy (Eq. [T]) is in general not exact and lifting it 
favors a "crystal" state in the energy landscape without 
modifying -if it remains sufficiently small- the dynamical 
aspects we have described. 

We have compared the results with the experiments 
on the kagome compounds. The present study gives a 
model for the spin freezing observed at Tg and provides 
an interpretation for the nature of the low-temperature 
phase. The picture of the "frozen" phase that emerges 
is that of an heterogeneous state with dynamical and 
frozen regions. The weak measured frozen moment is in- 
terpreted as the consequence of the remaining dynamics 
of the shortest loops and its strength is close to what 
is measured in the experiments. While in magnets in 
general, the on-site moment is reduced by the small os- 
cillations around the ordered state (spin- waves), here the 
main effect is argued to be the large-amplitude motion 
of the shortest loops. The short loop fluctuations do 
not fully destroy the moment for T < Tg but their pres- 
ence is in agreement with the persistent fluctuations ob- 
served by different experimental techniques (neutrons, 
/iSR, NMR). In particular, the observation in NMR of 
nuclei with different time-scales is consistent with the 
heterogeneous picture of the dynamics proposed here. 
In conventional magnets, the thermal excitations of the 
spin-waves destroy the on-site magnetization. Here, one 
needs longer loops, that are thermally excited only for 
T > Tg. These fluctuations give a spectral response that 
obeys a power-law ijj~^i^ in the small energy limit, very 
different from that of conventional magnets (flat response 
in two dimensions). A broad power-law response is in- 
deed observed experimentally in neutron inelastic scat- 
tering. Although the exponent seems to be underesti- 
mated, the experiments may not have had access to the 
low-energy limit or the exponent may be inacurately pre- 
dicted because of the interaction between the spin- waves 
and the discrete modes. In the paramagnetic phase, the 
model has algebraic spatial correlations at equilibrium 
(T > Tg), a feature that is not observed in elastic neutron 
scattering. We believe that this is not redhibitory, for the 
spin freezing we have described is not related to the long- 
distance behavior. In two spatial dimensions, the correla- 
tion length is always flnite at finite temperaturesiiii Fur- 
thermore, the chemical disorder is present to an amount 
which is difficult to quantify and which has been com- 
pletely neglected here. 

The energy scale that governs the freezing tempera- 
ture Tg is argued to be J in the small anisotropy limit 
(dynamically generated barriers), Tg = 0.04JiS' and it 
crosses over to Tg = 0.225DS^ in the strong anisotropy 
limit, typically if D/J > 0.18/5. This led us to a tenta- 
tive classification, where SCGO is in the small anisotropy 
limit and (Il30)Fc3 (804)2 (011)6 in the strong anisotropy 
limit. This is clearly a different interpretation from that 
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of chemical disorder, where Tg is governed by the amount 
of disorder <^ 

In order to disentangle intrinsic effects from the ef- 
fects of chemical disorder, one can test the present the- 
ory, in particular by characterizing experimentally the 
active magnetic degrees of freedom, for instance by neu- 
tron form factors)^ or by inferring the nanoscopic size of 
the frozen clusters. 
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